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ABSTRACT 


Theoretical argumer s for a 30% Increase In the solar luminosity over 
the past 4.7 billion years are reviewed. A scaling argument shows that this 
Increase can be predicted without detailed numerical calculations. The 
magnitude of the Increase Is Independent of nuclear reaction rates, as long 
as conversion of hydrogen to helium provides the basic energy source of the 
Sun. 

The effect of the solar luminosity Increase on the terrestrial climate 
is briefly considered. It appears unlikely that an enhanced greenhouse effect, 
due to reduced gases (NH-, CH, ) , can account for the long-term paleocllmatic 
trends . J 4 

INTRODUCTION 


Climatically significant changes of the solar luminosity (L) have been 
postulated to occur on time scales ranging from a few years to billions of 
years. The shorter time scales have been discussed extensively at this 
conference. In the present review, I will restrict myself to the longest time 
scales (S 10^ yr.) and discuss the basis for the astrophyslcal conclusion that 
the Sun was ~ 30% fainter 4.7 x 10’ yr. ago and that the evolution since the 
Sun's formation requires a slow, but steady, increase in L. 

I should note that this Is the only change in L predicted by stellar 
evolution theory, in its standard form. This prediction is common to all 
modern calculations and is supported by a large body of data from observational 
stellar astronomy (see reference 1 for a review of the observational evidence) . 
Nevertheless, the validity of this result has been questioned because of the 
apparent conflict with proxy indicators of the Earth's past climate (ref. 2-4). 
For this reason, a review of the theoretical arguments for the long-term 
increase of L is in order. 

Stellar evolution is governed by nonlinear differential equations derived 
from conservation laws and considerations of energy transport processes. 
Analytic solutions do not exist for any cases relevant to the Sun so numerical 
solutions must be used. Modern calculations require complex computer codes 
incorporating a variety of physical data on nuclear parameters, transport 
coefficients, and thermodynamic properties. In this respect, the situation is 
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similar to that encountered in current theoretical investigations of the 
terrestrial climate. This may, in fact, explain the reluctance of the climato- 
logists to accept the astrophysical i »ult: climatologists understand the 
pitfalls of accepting solutions obtained from complex computer codes at face 
value. For this reason, I will largely avoid discussion of numerical models 
and base the astrophysical case on simple t.aling laws. 


A SCALING MODEL OF THE SUN 


We begin by requiring that the Sun be in hydrostatic equilibrium, with 
gravitational forces balanced by the pressure gradient. The free-fall time 
of the Sun is on the order of an hour and any departures from hydrostatic 
equilibrium would show up as luminosity and radius changes on this time scale. 
For the spherically symmetric case, hydrostatic equilibrium is expressed as 


dP 

dr 



P 


( 1 ) 


where G is the gravitational constant, P and p are the pressure and density at 
a distance r from the center and m is the mass Interior to r. Measurements 
of the visible solar disk show that the Sun is spherical to within 1 part in 
10 5 (ref. 5). 

We can construct a one-zone model by replacing (1) by a finite-difference 
equation evaluated between the center c and surface s, with mean values 
enclosed in brackets ( ) : 
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Applying the boundary conditions P - 0, r « R (radius), and r - 0 gives 
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The scaling laws for the mean values are: 
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where m is the total mass. Inserting these scaling lavs into equation (2) 
gives 


P 

c 



( 6 ) 


To proceed further*, we necu an equation of state, relating the pressure 
to the density and temperature (T>. For typical conditions characterizing 
the bulk of the solar interior (;■ at 1 g/cm' and f at 10” to 10^ °K) , Coulomb 
interaction energies are at least 2 orders-of-magnltude smaller than particle 
kinetic energies. Thus, the ideal gas lav is an excellent approximation and 
this is what differentiates a star from a planet. Applying the ideal gas lav 
to the center gives 


P 


c 



(7) 


where k is the Boltzmann constant and the particle density is expressed is t».* 
Avogadro number h. divided by the mean mass w (in atomic mass units) par free 
particle. Eliminating P between (6) and ;7) and noting that the central 
density must scale as th$ mean density gives 


T * 
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If the scaling lav (5) is used to replace R, we get 
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Ue now turn to the question of how energy is transported from the core, 
where nuclear reactions produce energy, to the surface. Due to the high 
temperatures, radiative transport of energy is very efficient and dominates 
over the bulk of the interior. The mean free path i of a photon is typically 
1 cm so the photon diffusion approximation is valid to ordar t/R at 10“ l . 

The radiative diffusion equation viw.. spherical symMtry is 

i . 84ao r^T 3 dT /<•>• 
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where L is the total flux across a spherical surface at distance r from the 
center, r o is the Stefan-Boltzmann constant and k Is the Rosseland-mean opacity 
coefficient. Again, we use the one-zone difference approximation to write 
this equation as 



2_3 T - T 

■ .■irq , r T \ _c s . 

3 ' iep ' r - r 

c s 


Application of our previous boundary conditions plus T ■ 0 (i.e. T << T ) 

. S 8 C 

gives 


2 3 4 

. R T T RT 

L <x — 

Kp R icp 


( 10 ) 


where, since we are now dealing with a scaling law (proportionality), L 
can be replaced by L, r by R, etc. 

To evaluate the opacity coefficient we note that, in the solar 
interior, hydrogen and helium will be completely ionized and the heavier ions 
will be stripped of most of their electrons. Hydrogen and helium affect ic 
through free-free transitions while the heavier elements contribute primarily 
through bound-free transitions. Both processes are reasonably represented by 
the hydrogenic approximation so the absorption coefficient for a given ion 
varies inversely with the cube of the frequency. Although individual 
ionization states may contribute "noise" to the detailed dependence of ic on p 
and T, the broad dependence is given by Kramers’ oparity: 
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Putting this result into equation (10), and using equation (8) to eliminate 
the temperature, gives 


5.33 0.17 7.5 

m p v 


( 12 ) 


- 14 ^ e present rate of mass loss, due to the solar wind, is roughly 
10 mj yr. (ref. 6) and there is no reason to believe that the mass loss rate 

in the past was great enough to significantly affect m. The density 
dependence in (12) is so weak that we may also neglect changes in this 
parameter. Thus, the luminosity is primarily dependent on the mean molecular 
weight u and we rewrite (12) as 

L(t) - L(o) l^j-J (13) 
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We let X, Y, end Z denote the fractional abundances, by mass, of hydrogen, 
helium, and metala (X+Y+Z ■ 1). For a fully ionised gas, 


P « 


3 1 

2X + ft + yZ 


1 + 


3X + yY 


(14) 


The mean molecular weight Increases with time due to conversion of hydrogen 
(p “ 1/2) into helium (p ■ 4/3), by nuclear reactions, producing Q ■ 6x10*® 
erg per gram of hydrogen consumed. This energy must supply the luminosity 
of the Sun. Since Xm is the total mass of hydrogen, 


. w dX . _ L 
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Differentiating (14) with respect to time, and noting that 
dY/dt - - dX/dt, we get 

lE-.i 2 dX^ 5 m 2 l 
dt 4 p dt 4 mQ * 


(16) 


Finally, we can eliminate p between equations (13) and (16) . The resulting 
differential equation can be directly Integrated to give 


l <0 - MO) /l - ? t J - 15 '” 


(17) 


Since nuclear reactions are confined to the core, the present photo- 
spheric abundances should reflect the initial composition. Thus, we may 
evaluate p using X » 0.71 and Z ■ 0.02. Equation (17) becomes 


L 




1 -0.35 L(o) 


(tr m 


08) 


where L is expressed in units of the present solar luminosity v taken as 
3.9xlO J3 erg/s) and t is the present solar age (4.7x10* yr.). The initial 
luminosity required to match the present solar luminosity at t is L(o) - 0.76. 
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The scaling argumants predict that tha Sun was Initially 24X fainter than 
the present luminosity. A comparison of the luminosity evolution according 
to equation (18) with results from detailed numerical models (ref. 7) is 
shown in figure 1. As noted by 0. 0. Gough, the evolution predicted by 
numerical models is accurately represented by 


L ■ 



(19) 


where L is the present solar luminosity. This formula, rather than equation 
(18), it reconmended for studies of the evolution of planetary atmospheres. 


SUMMARY OF THE ASTROPHYSICAL CASE 


The above analysis shows that a quantitative prediction of the 
evolutionary increase of the Sun's luminosity may be made without detailed 
knowledge of the physical processes taking place in the interior. Therefore, 
this prediction is not affected by the uncertainties in this knowledge. In 
particular, we did not have to specify any nuclear reaction rates since the 
net reaction rate, Integrated over the solar mass, is determined by the 
measured solar luminosity. This is quite different from the case of the solar 
neutrino prediction, which is very sensitive to detailed nuclear reaction 
rates (ref. 8) . The discrepanc y between the predicted and observed neutrino 
flux should not be used to argue that the luminosity prediction is also 
questionable . 


APPLICATION TO THE EARTH'S CLIMATE 


Sagan and Mullen (ref. 9) pointed out that an enhanced greenhouse effect, 
due to higher concentrations of NH. and CH. in the Earth's atmosphere, could 
have maintained a warm climate even with a^lower solar luminosity. A similar 
conclusion was reached by Hart (ref. 10). This mechanism cannot, however, 
compensate for all of the solar luminosity evolution. 

Paleological evidence (ref. 11) shows that the Earth's atmospheric 
chemistry changed from reducing to oxidising some 1.5 to 2 billion years ago 
and this would have removed the enhanced greenhouse effect due to reduced 
compounds. Roughly one-half of the solar luminosity Increase occurs during 
the last 2 billion years but there is no evidence for a parallel Increase in 
the Earth's mean surface temperature. Indeed, isotopic studies of Precambrlan 
samples by Knauth and Epstein (ref. 12) indicate that the mean surface 
temperature has been decreasing during this time. Clearly, there is a need 
for further studies of the effects of crustal movements and volcanism. 
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biological activity, etc. on the long-term evolution of the Earth's climate. 
At present, it appears that the effects of solar evolution are still buried 
in the "noise" due to other uncertainties in paleoclimatic models. 
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Figure 1. Long-term evolution of the solar evolution. The 
evolution predicted by the scaling model (Equation (18)) is 
shown by the dashed line and the prediction from a detailed 
computer model is shown by the solid line. 




